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ABSTRACT 

The relative age of galaxies at different redshifts can be used to infer the Hubble parameter and put 
constraints on cosmological models. We select luminous red galaxies (LRGs) from the SDSS DR7 and 
then cross-match it with the MPA/JHU catalogue of galaxies to obtain a large sample of quiescent 
LRGs at redshift z ~ 0.03 — 0.39. The total 23,883 quiescent LRGs are divided into four sub-samples 
according to their velocity dispersions and each sub-sample is further divided into 12 redshift bins. 
The spectra of the LRGs in each redshift and velocity bin are co-added in order to obtain a combined 
spectrum with relatively high S/N. Adopting the GalexEV/StcLib model, we estimate the mean ages 
of the LRGs from these combined spectra by the full-spectrum fitting method. We check the reliability 
of the estimated age by using Monte-Carlo simulations and find that the estimates are robust and 
reliable. Assuming that the LRGs in each sub-sample and each redshift bin were on average formed at 
the same time, the Hubble parameter at the present time Hq is estimated from the age-rcdshift relation 
obtained for each sub-sample, which is compatible with the Hq value measured by other methods. We 
demonstrate that a systematic bias (up to ~ 20%) may be introduced to the Hq estimation because 
of recent star formation in the LRGs due to the later major mergers at z < 0.4, but this bias may 
be negligible for those sub-samples with large velocity dispersions. Using the age-redshift relations 
obtained from the sub-sample with the largest velocity dispersion or the two sub-samples with high 
velocity dispersions, we find Hq = 65^lkm s~ 4 Mpc" 1 or Ho = 74+4km s" 1 Mpc" 1 by assuming a 
spatially flat ACDM cosmology. With upcoming surveys, such as the Baryon Oscillation Spectroscopic 
Survey (BOSS), even larger samples of quiescent massive LRGs may be obtained, and thus the Hubble 
parameter can be measured with high accuracy through the age-redshift relation. 
Subject headings: cosmological parameters - cosmology: theory - galaxies:evolution - galax- 
ies: abundances - galaxies: stellar content 



1. INTRODUCTION 

The expansion history of the universe are presently 
studied with a few observational probes, such as the 
supernova la, baryon acoustic oscillations (BAO), weak 
gravitational lensing, and galaxy clusters, etc. Each of 
these probes has its pros and cons, and suffer from differ- 
ent s ystematic uncertainties (e.g., iFreedman fc Madord 
l2"010h . A new observational probe of the cosmic expan- 
sion history would be invaluable, and can provide ad- 
ditional cross check with the results obtained from the 
existing methods. Combining the results obtained by dif- 
ferent means may further help to constrain robustly the 
d ynamical nature of the universe . 

Uimenez fc Loebl (|2002| ) proposed a novel approach to 
explore the expansion history of the universe, which is 
based on the age-redshift relation of passively evolving 
massive galaxies. Assuming that the passively evolving 
galaxies at different redshifts were born approximately at 
the same time, the age of these galaxies can then be taken 
as a cosmic chronometer. If the ages of such galaxies can 
be accurately estimated, then this age-redshift relation 
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may be used to determine the cosmic expansion history. 
Even if there is some systematic errors in the absolute 
age measurements, it is argued that such errors could be 
canceled in the relative age of these galaxies at different 
redshifts, thus providing a good measurement of H(z): 



Indeed, observations show that the most massive 
galaxies are mainly composed of old stellar pop- 
ulations formed at redshifts z > 1 — 2, less 
than 1% of their present ste ll ar mass is formed 
at z < 1 (iDunlop et al.l 119961: iSpinrad et alj 119971; 
Cowie et all 119991: iHeavens et al.l l2004t IThomas et al.l 
20051: ICimatti et all 120081: iThomas etaLl 12010ft . hence 
t hese galaxies are suita ble for this application. 

Uimenez et al.l ((2003ft applied this method to a sample 
of massive galaxies at low redshift by fitting their spectra 
with the single stella r population (SSP) s pectra based on 
the SPEED model ([Jimenez et al. 2004), and obtained 
H = 69 ± 12km s" 1 Mpc" 1 . iSimon et all (|2005[ ) assem- 
bled a high redshift data set obtained from the Gemini 
Deep Survey (GDDS) and some other archival data, and 
applied the same method to estimate H(z) for a large 
redshift range (z ~ 0.1 — 1.8). These earlier works on the 
age-redshift relation adopted the SSP to fit each galaxy 
spectrum in the sample, and selected the age of the old- 
est one in each redshift bin as the envelop of the age. 
However, the poor signal-to-noise ratio (S/N) spectra of 
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individual galaxies and the contamination from the tel- 
luric emission and absorption may lead to uncertainties 
in the age estimates, and the method of the oldest galaxy 
envelop draw results from a small number of galaxies at 
the extremes of the distribution, which may also under- 
mine the validity of the result, and makes the method 
hard to use. 

To overcome this problem, I Carson fc Nicholl (|2010fi 
obtained the combined spectra for those luminous red 
galaxies (LRGs) with similar physical properties in each 
redshift bin by co-adding their spectra, of which the S/N 
is much higher than individual galaxies. They then esti- 
mated the age of the combined spectra by using the stan- 
dard Lick absorption line indices, which may be regarded 
as the mean age of a large sample of galaxies. They ob- 
tained the age-redshift relation, but they did not use this 
relation to further constrain the Hubble parameter. 

In this paper, we first select a LRG sample from the 
SDSS data release 7 (DR7). In order to improve the 
S/N and remove the contamination, we also use the 
combined spectrum rather than the single spectrum of 
each galaxy. However, we adopt the full spectrum fitting 
method, different f rom the standard Lick ab sorption line 
indices adopted by I Carson fe Nicholl (|2010t) . to estimate 
the mean age of the combined spectrum, and then obtain 
the age-redshift relation. Furthermore, we also use the 
age-redshift relation obtained from the combined spectra 
to constrain the Hubble parameter at the present time Hq 
and analyze the possible systematic bias in the estimated 
Hq. The paper is organized as follows. In Section 2, we 
describe the selection criteria of the LRG sample. In Sec- 
tion 3, we provide the details of the fitting method and 
the age-redshift relation estimated from the LRG sam- 
ple. In Section 4, we constrain the Hubble parameter by 
using the obtained age-redshift relation. Discussions on 
the resulted age-redshift relation and the possible associ- 
ated systematic bias are given in Section 5. Conclusions 
are summarized in Section 6. 

2. SAMPLE SELECTION 

In order to obtain the age-redshift relation and use it 
to measure the Hubble parameter, it is necessary to first 
select a large sample of passively evolving galaxies that 
contains the oldest populations with homogeneous phys- 
ical properties. The SDSS is currently the largest survey 
that provides hundreds of millions of detected objects 
with accurate photometric and astrometric c alibrations, 
and p art of the o b jects have excellent spectra (iPicr et all 
12003. iHogg et al.ll200l1 ). It is generally accepted that the 
LRGs are passively evolving galaxies and that they host 
the oldest stellar populations. Therefore, we pick our 
sample from the LRGs of t he SDSS DR7 (jYork et al.l 
l200q lAbazaiian et al. Il200l . 

For our purpose, it is necessary to determine the phys- 
ical properties of the LRGs with relatively high accuracy 
by using their spectra. Considering that some physical 
parameters such as the velocity dispersions and emission 
lines are available in the MPA /JHU sample only for CUT 
I LRGs (see lEisenstein et al.ll200l . we select only from 
the CUT I LRGs. To obtain accurate estimates of the age 
through the full spectrum fitting method, the selected 
galaxies should also have sufficiently high S/N. Our se- 
lection criteria for LRGs are similar to that described in 
I Carson fc Nichol (|2010l ). but with an additional restric- 
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Fig. 1. — The redshift distribution of LRGs. The solid histogram 
represents the distribution of LRGs in our final sample and the 
dotted histogram represents the distribution of all LRGs in our 
paper whose the lower-redshift peak is mostly contributed by the 
spiral bulges. 

tion on the S/N as follows: 

• selecting galaxies from Catalog Archive Server 
(CAS) database using the TARGET_GALAXY_RED flag; 

• selecting galaxies with S/N > 10 per pixel (in the 
continuum of the r-band wavelength range); 

• selecting galaxies which further satisfying the 
restrictions: specClass EQ 'SPEC_GALAXY',zStat 
EQ 'XCORRJiIC, zWarning EQ , eClass < 0, z 
< 0.4 and fracDevjr > 0.8 (see more details in 
I Carson fc Nicholl [2010h . 

According to the criteria above, 71,971 LRGs are se- 
lected from the SDSS DR7 and their redshift distribu- 
tion is shown in Figure Q] (the dotted histogram). These 
LRGs are probably contaminated by the bulges of late- 
type galaxies at low redshift, due to the limited fiber size 
(3") of the SDSS spectrograph. Such late contaminants 
would however have new star formations, and it is well 
known that the [On] and H a lines are indicators of star 
formation, hence we can remove the spiral bulges from 
the selected sample and obtain an sample of quiescent 
LRGs by using the spectral line data. Here we use the 
MPA/ JHU spectral line dat£0, which has been widely 
used in selecting quiescent galaxies in the literature. We 
select those LRGs as quiescent only if their H a and [On] 
line emission are consistent with zero at 2a level. With 
this criterion, we obtain 27,208 quiescent LRGs. We note 
here that using other emission lines to select quiescent 
ga laxies may result i n a sim ilar quiescent LRG sample, 
as I Carson fc Nichol (|2010T l pointed out that all other 
emission lines (such as the joint constraint of Hp and 
[Oin] or Nil and Sn) show similar zero-emission line dis- 
tributions. 

To estimate the Hubble expansion rate H(z) by using 
the age-redshift relation, it is important to select sam- 
ples of galaxies wi th homogeneous physi cal properties as 
demonstrated by iCrawford et al.l (|2010l ). For this rea- 
son, the quiescent LRGs selected above are divided into 
four sub-samples according to their velocity dispersions, 

5 http://www.mpa-garching.mpg.de/SDSS/DR7/raw_data.html 



The age-rcdshift relation For LRGs 



3 



TABLE 1 

The total number of galaxies in each sub-sample. 



Sample 




Velocity Dispersion Range 




Number 


sub-sampk 


i I 


200km s" 1 < ov < 230km s~ 


i 


4756 


sub-sample 


11 


230km s" 1 < <r v < 260km s" 


1 


8748 


sub-sample 


111 


260km s" 1 < ov < 290km s~ 


i 


7149 


sub-sample 


IV 


290km s" 1 < ov < 320km s~ 


i 


3230 


total 




200km s" 1 < o v < 320km s~ 


i 


23, 883 



which are listed in the MPA/JHU galaxy catalog. The 
velocity dispersion bins for each sub-sample of the LRGs 
are 200km s" 1 - 230km s -1 , 230km s" 1 - 260km s" 1 , 
260km s" 1 - 290km s" 1 and 290km s" 1 - 320km s" 1 , 
respectively. We denote these four sub-samples as sub- 
sample I, sub-sample II, sub-sample III, and sub-sample 
IV, respectively. In each of the velocity dispersion bins, 
the number of galaxies is still sufficiently large for the fol- 
lowing co-adding spectra to reach a high S/N (> 40). We 
do not consider galaxies with velocity dispersions larger 
than 320km s — 1 , as the total number of those galaxies at 
z < 0.14 is too small (< 30). In addition since the to- 
tal number of galaxies with velocity dispersion less than 
200km s _1 is 1843, much smaller than the number of 
galaxies in all other sub samples, we also exclude these 
galaxies. The number of galaxies in each sub-sample is 
listed in Table [1] and after excluding those LRGs with 
<j v > 320km s _1 or o~ v < 200km s _1 , the total number 
of the final sample is 23,883. The redshift distribution of 
these galaxies is shown in Figure Q] (the solid histogram). 

3. SPECTRAL FITTING METHODS 

There are three commonly adopted methods for mea- 
suring the age and mctallicity of a stellar system from its 
spectrum: (1) the SED fitting, (2) the Lick indices fit- 
ting, and (3) the full spectrum fitting. The first method 
is only sensitive to the general shape of the continuum, 
the second one focuses on using the strength or equiva- 
lent width of lines and specific spectrum features, and the 
third one accounts all the information of the spectrum, 
including both the continuum and the lines and specific 
features. The full spectrum fitting method has several 
advantages, such as being insensitive to extinction or flux 
calibration errors, and it is also not limited by the phys- 
ical broadening of lines since the internal kinematics is 
deter mined simultaneou sly with the population parame- 
ters dKoleva et al.l [20081) . though it is insensitive to the 
element ratio effects because of yet no available models 
about these. On the other hand, it is more sensitive to 
the wavelength range of the spectrum adopted in the fit- 
ting and the resolution of the spectrum compared with 
the fitting with the Lick indices. 

3.1. The ULySS software 

ULySS is an open-source software package developed 
by a group in Universitc dc Lyon, which implements the 
full-spectrum fitting to study physical properties of stel- 
lar populations. In ULySS, an observed spectrum is fitted 
by a model spectrum, adopting a linear combination of 
non-linear components, optionally convolved with a linc- 
of-sight velocity distribution (LOSVD) and multiplied 
by a polynomial function. The multiplicative polyno- 
mial is adopted to absorb errors of the flux calibration, 
Galactic extinction and other factors which may affect 



the shape of the spectrum. It minimizes x 2 value by 
the MPFIT function when matching an observed spec- 
trum with the model ones. The line spread function 
(LSF), an analogy to the point spread function (PSF) 
for images, is also introduced in ULySS in order to ef- 
fectively match the resolution of the model spectrum to 
the observatio ns. For detai l s abou t ULySS, the readers 
are referred to lKoleva etaLl (|2009al) . Since the full spec- 
trum fitting method may be sensitive to the wavelength 
range of the spectrum adopted in the fitting, we adopt 
the GalcxEV/StcLib model, a popular single stellar pop- 
ulations (SSPs) synthesis model which covers the largest 
wavelength range, i.e. 3200A-9500A. This wavelength 
range includes the Can triplet (AA8498, 8542, 8662A), 
which is a prominent feature produced primarily by an 
old population of red-giants and thus import ant for de- 
term ining the age of the quiescent LRGs (|Diaz et al.l 
[19891) . 

The GalexEV/SteLib(GS) population model is pro- 
duced by the isochron e synthesis code of BC03 
(jBruzual fc Charlotl[2003] ) . which is widely used in SDSS 
data analysis. It use the SteLib library which contains 
249 spectra, but only 187 stars have measured metal- 
licity and can be used to compute the predicted spec- 
tra with a 3Aspectral resol ution. The GS mod el adopts 
the Padova 94 i sochrones (iBertelli et al.l 119941) and the 
Chabrier IMF ([Chabrier fc Gillesl 120031) . Totally 696 
SSPs, covering the age of 0.1 ~ 20 Gyr and the [Fe/H] 
of —2.3 ~ 0.4 dex, are included in the GS model. The 
relevant information of the GS model is given in Table [2j 

3.2. Model pre-treatment 

The resolution match is a key issue in the model fitting 
because the resolution of a model spectrum is usually dif- 
ferent from that of the observational data. It is necessary 
to transform either the model spectrum or the observed 
spectrum to match the resolution of the other one. The 
spectral resolution is characterized by the instrumental 
broadening or the LSF. In practice, the LSF is not nec- 
essarily a Gaussian and m ay vary with wavelength (s ee 
more details on the LSF in lKoleva etafl (|2008l l2009al )). 
In ULySS, three types of calibrations (arc lamp, standard 
star, twilight spectrum) can be used to determine the 
relative LSF between the model and the observation. In 
this paper, we use the standard star to do the calibra- 
tions. ULySS adopts a convolution of the model with 
a series of LSFs, which arc determined at some wave- 
lengths, and then interpolates linearly in the wavelength 
range between two convolved models^ 

We follow the step s described in iKoleva et al.l (|2009a[) 
and lDu et al.l (|2010D to match the resolution of the spec- 
tra for galaxies in our sample with that from the GS 
model. We use the spectrum of SDSS standard star, 
which is already contained in the ULySS package, as a 
representative of the SDSS observations. For the GS 
model spectrum, the relative LSF is obtained in ULySS by 
comparing the spectrum of the observed spectrum of the 
SDSS standard stars with the GS model spectrum of the 
stars with the same physical properties. We then adopt 
this relative LSF to generate the resolution-matched GS 
model spectrum by the LSF convolution function in the 
ULySS package. Below when we mention the GS model 
spectrum we are actually referring to such resolution- 
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TABLE 2 
The settings in the GS model. 



Model Library Resolution Wavelength Age Z IMF Track 
(A) (A) (Gyr) (dex) 

GS SteLib 3 3200-9500 0.1-20 -2.3-0.4 Chabrier Padova 94 



matched ones. 

3.3. Spectrum fitting 

We adopt the SSPs given by ULySS to fit the combined 
spectrum or the spectrum of each galaxy in our sample 
but do not consider the detailed star formation history 
(SFH) of each galaxy, as the sample is almost homoge- 
neous and the galaxies in it are passively evolving. We 
defer the discussion of the effect due to the difference in 
the SFH of each galaxy to Section [5j 

3.3.1. Combined Spectrum 

High S/N spectra are essential to obtain accurate es- 
timates of the age of galaxies. We have tested the effect 
of the S/N on the age uncertainty. We found that the 
uncertainties in age estimates are about 20% and 10% for 
galaxy spectra with S/N = 20 and 40, respectively. In 
order to apply the age-redshift relation effectively, the 
uncertainties in age estimates need to be smaller than 
10%. For the majority of galaxies in our sample, how- 
ever, their SDSS spectra have S/N < 20, which leads 
to > 20% uncertainties in age estimates. To overcome 
this problem, we choose to co-add the spectra of galaxies 
in each sub-sample (and each red-shift bin) to obtain a 
combined spectrum with significantly higher S/N (> 40). 

To do this, we first correct the foregrou nd Galactic ex- 
tinctio n by using the reddening maps of Schlcgcl ct al. 
(|1998f ) for each galaxy. As the number of galaxies at 
z < 0.03 is too small in all of the sub-samples, these 
galaxies are neglected in the following analysis. The 
galaxies in each sub-sample are then divided into 12 red- 
shift bins from z = 0.03 to z = 0.39 with a redshift step 
of Sz = 0.03. A combined spectrum is then obtained for 
each redshift and velocity dispersion bin by co-adding the 
rest frame spectra of all the galaxies in that bin through 
the IRAF task SCDMBINE. We thus obtained 12 combined 
spectra over redshift bins from 0.03 to 0.39 for each sub- 
sample. 

The combined spectra are then fitted by ULySS. ULySS 
use Levenberg-Marquardt routine to search the param- 
eter space to get the minimization of x 2 - So it needs 
some initial guess value to begin searching the parame- 
ter space. In this paper, we set the initial guess for the 
age to 8 Gyr and for metallicity to solar metallicity since 
the LRGs are expected to be old and metal rich. We 
do not set any limit on the allowed age and the abun- 
dance of metallicities of the model spectra. The fitting 
is performed in the whole wavelength range covered by 
the combined spectra. The estimates on the age, metal- 
licity and velocity dispersion of the model galaxies and 
the associated errors are then obtained for the best-fit 
model to the combined spectrum. For illustration, Fig- 
ure [2] shows the combined spectrum in the first redshift 
bin of the sub-sample II, and the best-fit model spectrum 
by the GS model and its residues. The fitting results for 
the sub-sample II are listed in Table [3] 



Figure [3] illustrates the best fit result obtained by 
adopting the GS model. The top, middle and bottom 
panel in the Figure shows the velocity dispersions, metal- 
licities and ages for each redshift bin in each of the four 
sub-samples, respectively As seen from the top panel of 
Figure El the fitting results on the velocity dispersions 
are consistent with the ones given in the MPA/JHU cat- 
alogue, which are shown in Figure [3] as open squares, 
triangles, diamonds, and circles, respectively. The mid- 
dle panel shows that the galaxies in each sub-sample 
have similar metallicities, confirming that our samples 
are almost homogeneous, though of course subtle differ- 
ences remain. From the bottom panel, it is clear that 
the mean age (i age ) decreases with rcdshifts from 0.03 
to 0.39. Apparently, the age of the galaxies with higher 
velocity dispersion tend to be somewhat older than those 
with lower velocity dispersion, which is consistent with 
the well known "downsizing" formation of galaxies, i.e., 
the bigger and more massive galaxies were formed ear- 
lier, while the small galaxies formed later (jCowie et al.1 
119961) . Assuming that the < agc — a v relation follows a 
power- law, i.e., < ag0 oc crj, we fit the relationship for each 
redshift bin and obtain the slope 7. The mean value of 
7 for all the redshift bins (and its standard deviation) 
is ~ 0.77 ± 0.25. We may also first obtain the mean 
iage for each sub-sample and then fit the mean i agc — a v 
relationship and find 7 ~ 0.79 ± 0.24. 

A number of studies have obtained the relationship 
between the age (i age ) and velocity dispersion (a v ) of 
early- type galaxies or L RGs and shown t hat th e t age de- 
pends a v . For example, iCaldwell et al.l (|2003| ) analyzed 
the integrated spectra of 175 nearby early- type galaxies 
by using higher order Balmer lines as the age indicators 
and found that early-type galaxies with lower <r v have 
smaller lu minosity-weigh t ed me an ages. The data ob- 
tained by ICaldwell et al.l (|2003| ) suggests th at the slope 
of the £ are — cr„ relation is ro ughly 0.8—1.2 (|Nelan et al.l 
I2005h . IThomas et "all (|2005| ) studied the spectra of 124 
early-type galaxies in both high and low density environ- 
ments by using the absorption line indices, and they also 
found that t ap - R corre lates with a v . Adopting the data in 
IThomas et al.l (|2005| ) , th e slope of the t nne . — a v relation is 
found to be 0.78 ± 0.23 (jNelan et al.ll2005h . INelan et al.l 
(|2005h investigated the spectroscopic line strength of 
4097 red-sequence galaxies in 93 low-redshift galaxy clus- 
ters, and they found i agc oc (j0- 59 ± - 13 . iSmith et al.1 
(|2009t ) studied the spectra of 232 quiescent galaxies in 
the Shapl ey supercluster and found t agc oc cr°- 40 . As dis- 
cussed in INelan et alj (|2005h , the differences in sample 
selection criteria and emission treatment may be account 
for the differences of i agc — c scaling relation. Consid- 
ering of these differences, our results are well consistent 
with those obtained by previous works. 

3.3.2. Reliability of the fitting results 
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Fig. 2. — The combined spectrum of the galaxies in the first redshift bin of the sub-sample II and its best-fit by adopting the GS model. 
In the top panel, black and blue lines show the combined spectrum and its best-fit model spectrum, respectively. These two spectra are 
almost superimposed on each other and the black line can be seen only when zooming in the Figure. The cyan line is the multiplicative 
polynomial to absorb the effects of an imprecise flux calibration and of the Galactic extinction. The red regions are rejected from the 
fitting (rejection of flagged telluric lines and automatic rejection of outliers). The bottom panel shows the residuals of the best-fit. The 
solid green lines mark the lcr deviation, and the dashed green line represents zero residuals. 



TABLE 3 

Results of a SSP fitting with the GS model for sub-sample II. 



Redshift Interval 


SSP-equivalent age 


SSP-equivalent [Fe/H] 


velocity dispersion 




(Gyr) 


(dex) 


(km s" 1 ) 


0.03 < z < 0.06 


7.13±0.53 


0.17±0.01 


243.5±4.0 


0.06 < z < 0.09 


6.52±0.43 


0.17±0.01 


239.2±3.7 


0.09 < z < 0.12 


6.31±0.40 


0.16±0.01 


236.0±3.7 


0.12 < z < 0.15 


6.01±0.40 


0.17±0.01 


239.0±3.7 


0.15 < z < 0.18 


6.11±0.45 


0.15±0.01 


240.7±3.9 


0.18 < z < 0.21 


5.49±0.24 


0.18±0.01 


242.4±3.8 


0.21 < z < 0.24 


5.37±0.22 


0.17±0.01 


243.0±3.5 


0.24 < z < 0.27 


4.98±0.22 


0.15±0.01 


242.7±3.5 


0.27 < z < 0.30 


4.96±0.19 


0.14±0.01 


247.6±3.3 


0.30 < z < 0.33 


4.78±0.16 


0.12±0.01 


247.4±3.4 


0.33 < z < 0.36 


3.82±0.24 


0.16±0.02 


243.7±3.5 


0.36 < z < 0.39 


3.64±0.19 


0.17±0.02 


248.9±3.2 


Note. — Column 1 


is the redshift interval, columns 2, 3 and 4 are 


the SSP-equivalent age 



in unit of Gyr, the [Fe/H] in unit of dex, the velocity dispersion in unit of km s 1 and their 
associated errors, respectively. 



It is important to check the reliability of the fitting 
results to ensure it is not highly dependent on the par- 
ticular synthesis technique adopted, or affected by the 
existence of multiple solutions due to degeneracies among 
t he model parameters . 

iKoleva et al.l (|2008t ) analyzed the spectra of Galactic 
clusters using ULySS, and they found that stellar popula- 
tions of these clusters obtained from the model are well 
consistent with that ob t ained f rom the color-magnitude 
diagrams. IKoleva et al.l (|2009bD further analyzed the de- 
tailed star formation history of 16 dwarf galaxies by using 
either ULySS or STECKMAP, and they found that the two 
programs give remarka bly consistent results. In addition, 
iMichielsen et al.l ()2007l ) adopted two different techniques, 
i.e., the Lick/IDS index system and the full spectrum 
fitting method, to test the reliability of the estimates of 
the ages and metallicities of 16 dwarf elliptical galax- 
ies, and they found these two techniques give consistent 
results on the age and the metallicity, with an rms er- 



ror of 1.63 Gyr in age and 0.09 dex in [Z/H]. IDu et all 
(|2010f ) synthesized the star formation histories and evo- 
lution of 35 brightest E+A galaxies from the SDSS DR5, 
and demonstrated the robustness of the ULySS technique 
in measuring the age and metallicity of stellar systems. 
These studies show that the fitting results with different 
synthesis technique are consistent with each other, and 
the results of our ULySS fitting should be robust. 

We perform Monte-Carlo simulations to visualize the 
degeneracies and validate the errors by simulating the 
effect of the noise. For each combined spectrum, we per- 
form 1000 simulations. In each of the simulations a ran- 
dom Gaussian noise is added to the combined spectrum 
and the amplitude of the added noise is set to the esti- 
mated noise associated to the combined spectrum. We 
then get the mean values of the age, abundance of metal- 
licity and velocity dispersion, and their standard devia- 
tions, correspondingly. Figure [4] shows the results of the 
1000 Monte-Carlo simulations (open symbols) and the 
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Fig. 3. — The mean galactic properties extracted from the combined spectra by using the GS model. Panels (a), (b), (c) and (d) show the 
velocity dispersions, mctallicities, mean luminosity-weighted ages and mean formation time of the LRGs obtained through their combined 
spectra in each redshift bin for each of the four sub-samples, 200km s _1 < ov < 230km s~ 1 (purple squares), 230km s _1 < cr v < 260km s 
(red triangles), 260km s~ 1 < ov < 290km s^ 1 (green diamonds), and 290km s~ 1 < o v < 320km s — 1 (blue circles), respectively. For clarity, 
those symbols for the same redshift bin but different sub-sample are put offset by 8z = —0.002, 0, 0.002, and 0.004 along the horizontal 
axis for the four sub-samples with velocity dispersions from low to high, respectively. The galactic velocity dispersions obtained directly 
from the MPA/JHU catalogue are shown as open squares, triangles, diamonds, and circles for each sub-sample, respectively. In the panel 
(c), the dashed line indicates t\j(z) — 5.5 Gyr, where t\j(z) is the age of the universe for a ACDM cosmology and is for reference only. The 
concordant cosmological model, i.e., Ho = 71km s -1 Mpc — 1 and f2 m = 0.27, is adopted to obtained the mean galaxy formation time shown 
in panel (d). 



original best fit shown in Figure [3] (filled symbols). Ac- 
cording to Figure IH we conclude that in most of cases, 
the value of mean velocity desperation and mean age for 
every bin of the 1000 Monte-Carlo simulations lead to the 
best fitting solutions, though there are still some devia- 
tions especially for sub-sample IV, which may be caused 
by relatively small number of galaxies used to obtain the 
combined spectrum. Since we only need to model the 
age-redshift relation, the diversity of mctallicity between 
the best fitting value and Monte-Carlo simulations value 
do not affect our conclusion. 

The general agreement of the above series tests sug- 
gest that the ULySS technique is robust in determining 
the age, metallicity and velocity dispersion of stellar sys- 
tems, and the fitting results on the physical properties of 
galaxies obtained from the full spectrum fitting are also 
secure. 



4. CONSTRAINTS ON THE HUBBLE PARAMETER 

The age-redshift relation obtained from observations 
can be directly fitted by 



(2) 



where i age is the mean age obtained from the combined 
spectrum, if orm is the mean formation time of the qui- 
escent galaxies and assumed to be a constant for each 
sub-sample, and t\j is the age of the universe. According 
to the ACDM cosmology, tu at redshift z is given by 



1 

tv = w 



dz 



(3) 



which depends not only on the value of Hq but also the 
composition of the universe, i.e., Q m and 

In order to obtain a model-independent measurement 
of Hq, however, we first simply assume H(z) = Hq + 
H'z (but it may be a good approximation only at low 
redshift), then the age of the universe is given by 



tu = 



dz 



(l + z)(H + H'z) 



(4) 



Using the standard \ 2 minimization, we fit the age- 
redshift relation obtained from each sub-sample to 
get the best fit of Hq and fform, and the uncer- 
tainty of Hq by marginalizing over H' and tf orm 
(and the uncertainty of if or m by marginalizing over 
H' and Hq). Figure [5] shows the best fit of 



The age-rcdshift relation For LRGs 



7 



320 "^ 
■p 290 — 

CO 

£ 260 — 

b > 230 - 
200 _ 
— 0.25 E- 

X - 

D Z 

3 0.20 r 
^ 0.15 r 
^ 0.10 — 





0) 

o 



4 — 




6 — 

4 — 
0.00 




0.30 



0.35 



0.40 



Fig. 4. — Legends are similar to that of Figure [3] except that the Monte-Carlo simulation results (open symbols) are additionally 
shown here. And the Monte-Carlo simulations are performed to visualize the degeneracy and validate the estimation of the errors (see 
Section PSCZ) . 

the age-redshift relation for each sub-sample, respec- 
tively. The best fits of the Hubble parameter Hq 

Mpc" 1 , 77i 7 km s" 1 Mpc" 1 , 



range from 84*gkm s 

2+5l 



TABLE 4 

The best fit of the Hubble parameter at the 



68I 7 km s" 1 Mpc" 1 to 6314km s" 1 Mpc" 1 for the sub- 
samples with velocity dispersions from low to high (see 
Table 2]) , but the mean galaxy formation time can not 
be well constrained. The age-redshift relations obtained 
from all the four sub-samples can also be fitted simul- 
taneously, and only six parameters are now involved in 
the fitting, i.e., Hq, H', and the mean formation time of 
the galaxies in each sub-sample, denoted as iformi, £form2, 
£form3) and if rm4, respectively, since Hq and H' are the 
same for all the age-redshift relations. By marginalizing 
over parameters H', tf orm i, *form2, t lorm3 and t {oim4 , we 
obtain the best-fit of Hq = 73i 3 km s" 1 Mpc" 1 (see the 
left panels of Figure [6]) 

Assuming spatial flatness for simplicity, we further 
adopt equation [5] and equation [3] (given by the standard 
ACDM model) to fit the age-redshift relation obtained 
from each sub-sample, separately. The best fit of Hq 
is obtained by marginalizing over f2 m and <f orm and it 
ranges from 89tgkm s" 1 Mpc" 1 , 83tgkm s" 1 Mpc" 1 , 

Mpc -1 for the four 



to 65_2km s 



72I?km s" 1 Mpc~ 

sub-sample velocity dispersion from low to high, respec- 
tively (see Figure [5] and Table 0}. We also obtain Hq 
by fitting the age-relations obtained from all the four 
sub-sample simultaneously by marginalizing over J7 m and 
the mean formation time of those galaxies, and the best 
fit of Hq is SOi^km s _1 Mpc -1 (see the right panels of 



present time Hq 


Sample 


H(z) 


= H + H'z 


Flat ACDM 


H 


\ V 


Ho 




Sub-sample I 


84t£ 


1.47 


89±l 


1.43 


Sub-sample II 


77±? 


1.26 


83l 9 8 


1.21 


Sub-sample III 


68 + | 


0.72 


72 + ? 


0.66 


Sub-sample IV 


63±I 


0.86 


65±l 


0.86 


Sub-sample III+IV 


68+1 


0.82 


74±l 


0.78 


All sub-samples 


73+4 


1.14 


80^ 


1.08 



Note. — Here Hq is in unit of km s Mpc , x v ls the 
reduced \ . Columns 2 and 3 list the best-fit value of Hq, its 
lcr error and the reduced x 2 > by assuming H(z) — Hq + H f z; 
while columns 4 and 5 list the best-fit value of Ho, its la error 
and the reduced x 2 by assuming a spatially flat ACDM model. 

Figure [6]). And the best fit of f2 m is 0.09 and it ranges 
from 0.08 to 0.25, which cannot be constrained with high 
accuracy. By marginalizing over other parameters, we 
obtain the best fit of t foiml , tf or m2, £form3, or t fmmA as 



9.7 



+0.3 
-2.7 



Gyr, 



q 2+ 03 



Gyr, or 8.4+°J Gyr. 



Gyr, 8.8±§-J 

These values of tf orm seem to be far too large than 
the expectation from the standard ACDM cosmology 
for LRGs, which are mainly caused by the poorly con- 
strained fi m (=0.09). By assuming a flat universe with 
£l m =0.27 (the concordant cosmological model) and re- 
fitting the age-redshift relations for the four velocity 
dispersion bins simultaneously, the best fit gives Hq = 
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Fig. 5. — The best fits to the age-redshift relations obtained from each sub-sample (points with errorbars) by assuming H(z) = Ho + H'z 
(dashed lines) and a spatially flat ACDM model (solid lines), respectively. The dashed lines almost overlap the solid lines. 



73km s 1 Mpc l , t form i = 5.9 Gyr, i/„ m2 = 5.4 Gyr, 
tform3 = 5.0 Gyr, and tf ormi = 4.6 Gyr, respectively. 
Clearly the mean formation time of those less massive 
galaxies is smaller than that of those more massive galax- 
ies, which is fully consistent with the "downsizing" evo- 
lution nature of galaxy formation. One may also directly 
obtain the formation time from the age-redshift rela- 
tions by adopting the concordant cosmological model, 
i.e., fixing the Hubble constant Hq = 71km s _1 Mpc -1 , 
the matter density Q m = 0.27 and dark energy density 
Oa = 0.73, the formation time from the age-redshift 
relations are tformi = 6.2 Gyr, tf or m2 = 5.7 Gyr, 
£form3 = 5.3 Gyr, and tf orm 4 = 4.9 Gyr, respectively. 
The mean galaxy formation time is about 5.5 Gyr for 
the four velocity dispersion bins, which is adopted in Fig- 
ure [3] (represented by the dashed line). Note that the ob- 
tained tform is the average age of a population of galaxies 
but neither the age of individual galaxies nor the oldest 
population of stars in those galaxies. Note also there is 
strong degeneracy between tf orm and (iToi^m) obtained 
from the fitting, which introduces a large uncertainty in 
the estimation of tform- 

According to the above fittings, obviously a strong 
constraint on Hq can still be obtained by either assum- 
ing a simple model independent form of the evolution of 
H(z) or a flat universe, i.e., f2 m + Oa = 1, although the 
age-redshift relations are obtained in a limited rcdshift 
range and the uncertainties in the age estimates may 
be substantial. The Hq estimated from the sub-sample 
with lower velocity dispersion tends to be higher than 
that from the sub-sample with higher velocity disper- 



sion (more massive and luminous LRGs), which may be 
due to some bias introduced by the systematical differ- 
ence in the assem bly history of l ess ma ssive galaxies and 
massive galaxies. iBrown et al.l (|2007l ) pointed out that 
the evolution of galaxies in the red sequence is heavily 
dependent on luminosity, i.e., the lower the luminosity 
of the galaxies, the more significant the population of 
new stars formed since z = 1. So the contamination 
from the population of stars formed at low rcdshift (e.g., 
z < 0.4) is probably more significant in the sub-sample I 
and sub-sample II than that for the sub-sample III and 
sub-sample IV. And the age-redshift relation estimated 
from the lower velocity dispersion sub-sample tends to 
be shallower than that from the higher velocity disper- 
sion sub-sample, which may lead to an overestimation of 
the Hubble parameter up to ~ 20% (see discussions in 
Section I5.3|) . 

As the possible systematic bias may be not signifi- 
cant for the two subs-samples with high velocity dis- 
persions, we also fit the age-redshift relations obtained 
from the two sub-samples simultaneously by assuming 
either a spatially flat ACDM model or H(z) = H + H'z. 
The best fit of Hq is either 74tlkm s" 1 Mpc -1 or 

68tskm s" 1 Mpc" 1 . And the best fit of O m is 0.071^ 
if assuming a spatially flat ACDM model. 

Considerable progress has been made in determining 
the Hubble parameter over the pa st two decades by us- 
ing m any different techniques (e.g.. lFreedman fc Madorel 
120101 ). For example, the Hubble parameter is esti- 
mated to be Hq = 73 ± 5km s _1 Mpc -1 by using 



The age-redshift relation For LRGs 



0) 

cn 
o 




0.0 



0.2 

Z 



0.3 



0.0 



0.2 
Z 



0.3 



0.4 



100 
^ 90 
Q- 80 

T M 70 
E 

60 

o 
X 

50 



40 




I 00 

90 



a 80 



: to 

E 



30 40 



50 



60 70 
H' 



80 90 100 



70 
60 

50 

40 L 
0.0 




0.2 



0.4 



0.6 



Fig. 6. — The best fit to the age— redshift relations by fitting these relations obtained from the four sub-samples simultaneously. Left 
panels show the results obtained by assuming H(z) = Hq + H'z, while right panels show the results obtained by assuming a spatially flat 
ACDM model. The top panels show the best fits (lines) to the age-redshift relations (points with errorbars) obtained from the sub-sample 
I (purple), sub-sample II (red), sub-sample III (green) and sub-sample IV (blue), respectively. The left-bottom panel show the confidence 
levels of the fitting parameters (Hq,H') by marginalizing over the formation time tf oim i, tf orm 2, iform3i tform4i an d the right bottom panel 
shows the confidence levels of the fitting parameters (HojQm) by marginalizing over tf OIm i, tf or m2> *form3> *form4- The crosses denote the 
positions of the best fits, the contours correspond to the la, 2a, 3a levels, respectively. 

In this paper, we have obtained the mean age, mctal- 
licity and velocity dispersion for a sample of quiescent 
galaxies selected from the SDSS DR7 in different red- 
shift bins. The age-redshift relation derived from those 
quiescent galaxies is consistent with the expectation from 
the ACDM cosmology and may provide a good estimate 
of the Hubble parameter {Hq). However, the estimate 
of Hq is valid only if those quiescent galaxies are pas- 
sively evolving and the galaxies in different redshift bins 
represent the same population formed more or less at the 
same time. In order to check whether these requirements 
are satisfied for the quiescent galaxy sample selected in 
this paper, we perform some tests below to investigate 
the evolution effects due to different star formation his- 
tory or galaxy mergers, and illustrate that those quies- 
cent galaxies are indeed more or less formed at the same 
epoch using the evolution of their average colors. 



the tip of the red giant branch as an alternate cali- 
bration to the Cepheid distance scale (jMould fc Sakail 
I2008D ; 72 ± ^{random) ± 11 (systematic)km s" 1 Mpc" 1 
by using the surface b rightness fluctuation to deter- 
mine cosmic distances (jBlakeslee et al.l 120021) ; 74.3 ± 
3.6km s _1 Mpc -1 by using Type la supernovae 
(|Riess et al.l (20091 ); 76.9;^ts° n km s" 1 Mpc" 1 by usin g 
the Sunyaev-Zeldovich effect (jBonamente et ail I2006D ; 
and 68 ± 2km s _1 Mpc -1 by using the BAO sig- 
natur e in the matter power spectrum (jPercival et al.l 
l2010f >. The Hubble Space Telescope key project 
yielded a consistent value of Hq = 72 ± 3 (random) ± 
7(systematic)]im s _1 Mpc -1 by co mbining the data ob- 
taine d f rom different te c hniqu es (jFreedman &; Madord 
120101 ). iKomatsu et al.l (|2009l ) also obtained a value 
of Hq = 70.5 ± 1.3km s^ 1 Mpc^ 1 by combining the 
WMAP-5 data with t he SNe la and BAO data, while 
iTammann et al.l (|2008l ) found consistently low values of 
Hq, from several different tracers they obtain a mean 
value of Hq = 62.3 ± 1.3km s _1 Mpc -1 . Considering 
of the possible systematic error in the Hq estimation by 
using the age-redshift relation, our estimates of Hq are 
fully consistent with those listed above. 



5. DISCUSSION 



5.1. Star formation history 

According to the selection criteria, most galaxies in 
our sample should be quiescent galaxies and supposed to 
be passively evolving. However, the star formation his- 
tory of those galaxies may not be a single burst. To test 
whether there are significant younger stellar populations 
in the galaxies, we re-do the full spectrum fitting for each 
combined spectrum by adopting two stellar components, 
one young stellar population (YSP) and one old stellar 
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population (OSP). The age of the young stellar popula- 
tion is assumed to be in the range from 0.1 Gyr to 1 Gyr, 
while for the old population it is assumed to be in the 
range from 1 Gyr to 15 Gyr. For both populations, the 
metallicity is a free parameter without restrictions. Ac- 
cording to the fitting results, we note here that the fitting 
age of the YSP always reaches the edge of its limits (the 
age of YSP is always either 0.1 Gyr or 1 Gyr), which 
means that a YSP with an age in the given range can 
not be found. Furthermore, even if there exist a YSP, 
the light fraction (LF) is very small (2.8% on average) 
so that it can be neglected. Therefore, we conclude that 
the YSP (the age of which <1 Gyr) is negligible and not 
required in the fitting. We also test the cases by setting 
a larger upper limit on the age of the YSP, e.g., 2 Gyr, 
3 Gyr or even 4 Gyr (or alternatively a lower limit on the 
age of the old stellar population, i.e., 4 Gyr to 6 Gyr), 
and also find that no significant YSP is required by the 
fitting. All these tests suggest that most of the quies- 
cent LRGs may be passively evolving and not experience 
significant recent (< 2 Gyr) star formation. For those 
galaxies at low rcdshift bins, however, we find there may 
exist YSP's with age ~ 3 Gyr, possibly due to the later 
major mergers (see discussions in Section 15. 3|) . 

To close this sub-s ection, we note here that 
iToieiro fc Percivall (|2010t) find bright LRGs bein g consis- 
tent with pure passive evolution while faint LRGs slightly 
deviating from pure passive evolution as revealed by the 
evolution of the number and luminosity density of LRGs 
as well as that of their clustering. The lesser passive- 
ness of LRGs with smaller velocity dispersion seems not 
to be able to be directly revealed by the combined spec- 
tra of LRGs studied in this paper, which might be due to 
that the signature of YSPs (probably with quite different 
ages) in (some of) the LRGs may be diluted or smoothed 
due to the co-adding of a large number of LRG spectra 
in our analysis. High S/N spectra of individual LRGs 
may be helpful to clarify the less passiveness of small 
LRGs, however, the S/N of most LRGs in our sample 
are only slightly larger than 10 and not sufficiently high 
for clarifying this problem. 

5.2. Combined spectrum vs. single spectrum 

The combined spectrum in each redshift and velocity 
bin may only represent the mean spectrum of galaxies in 
that bin. Docs the physical properties derived from the 
model fitting of this combined spectrum represent the 
mean properties of all galaxies in that bin? In order to 
test this, we fit each single spectrum for all the galaxies in 
each velocity dispersion bin and redshift bin with ULySS. 
The initial settings are the same as that in Section IS. 3. II 
After obtaining the best fit of the age for each single spec- 
trum, we calculate the mean of the ages of those galaxies 
which belong to the same velocity dispersion and redshift 
bin and its standard deviation. Figure [7] shows the mean 
age-redshift relation for every single spectrum. As seen 
from Figure the mean age-redshift relation obtained 
from the single spectrum fitting is still consistent with 
that obtained from the combined spectra but with much 
larger errors. 

5.3. Major mergers of quiescent galaxies 

The differences among the model spectra become small 
if the ages of the stellar populations are larger than 



12 



10 - 



6 8- 
" 6- 



"ikf-J. fr-§Nl 



0.0 



0.1 



0.2 

z 



0.3 



0.4 



Fig. 7. — The mean age-redshift relation derived by fitting each 
single spectrum in the LRG sample with ULySS for each redshift 
bin of each of the four sub-samples, 200km s _1 < <r v < 230km s _1 
(purple squares), 230km s _1 < cr v < 260km s _1 (red trian- 
gles), 260km s -1 < <t v < 290km s -1 (green diamonds), and 
290km s _1 < ct v < 320km s" 1 (blue circles), respectively. The 
dashed line indicates % (z) — 5.5 Gyr , where t\j(z) is the age of the 
universe for a ACDM cosmology and is for reference only. 

4 Gyr. Therefore, the uncertainties in the estimates of 
the ages and metallicities of galaxies with low S/N spec- 
tra are substantial. For this reason, we have combined 
the spectra of galaxies in each redshift and velocity dis- 
persion bin together to improve the S/N. However, the 
age difference obtained by modeling the combined spec- 
tra can represent the age difference of the universe only 
if those galaxies at different redshifts are the same pop- 
ulation and were formed more or less in the same epoch. 
There is a potential caveat to this approach, i.e., if a 
significant fraction of quiescent galaxies experience ma- 
jor mergers and thus some star formation at low redshift 
z < 0.4, then some quiescent galaxies in the redshift bin 
0.03 — 0.06 were formed at z < 0.4 and they were not rep- 
resented by those galaxies in the redshift bin 0.36 — 0.39. 
It is important to check the effect of the major merger of 
galaxies on the age-redshift relation. 

The major merger rate of galaxies is defined as the 
number of mergers per galaxy with mass larger than a 
threshold (Af») per unit time and denoted as dN mrg / 'dt, 
and it c an be roughly e s timat ed by the fitting formula 
given bv lHopkins etaLl (|20Tot) . i.e., 



dN n 



at 



A(M min )(l + z)^ M ^[per galaxy] 



(5) 



where the normalization is 

A(M min ) major « 0.02[1 + (M min /M ) a5 ]Gyr- 1 , 
and the rcdshift evolution is 



/3(M min ) 



mm ) major 



1.65- 0.15 log(M min /M ), 



and Mo = 2 x 10 10 Mq. For the galaxi es in our sample, 
they are generally brighter than 3L* (|Eiscustci n et al.l 
120011) and have stellar masses in the range from lO n M0 
to a few times 1O 12 M according to iTal et all (|201lD . 
Assuming that the minimum mass of those galaxies is 
Mmm = 10 11 Mq, the average number of major mergers 
experienced by a galaxy at redshift z since z = 0.40 is 



NmrgM 



dNmrg | dt 

dt dz 



(6) 
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where dt/dz is given by Eq. ([T]). The fraction of galaxies 
in a rcdshift bin z±dz which experienced major mergers 
since z = 0.4 is - N mTg (z), and N mlg (z) ~ 0.33, 0.26, 
and 0.01 at z = 0.03, 0.11, and 0.39 for our sample, re- 
spectively. Note that almost all the major mergers are 
dry mergers for massive galaxies with mass > lO 11 -/^ 
similar to that in our sample (Hopkins et al. I I2010D . 
These dry mergers may also lead to new star formation 
or star burst in the galactic centers, and thus introduce 
a systematic bias to the age-redshift relation, because 
the mean age obtained from the combined spectra by 
the GS model at low redshift bins may be systematically 
underestimated due to the additional population of stars 
formed later. 

We check the effect of these major mergers on the age- 
redshift relation as following. First, we assume an under- 
lying age-redshift relation according to the ACDM, i.e., 
the age of the galaxy at redshift z is t\j(z) — if orm , and 
Hq = 71km s" 1 Mpc" 1 , f2 m = 0.27, i form = 5.5 Gyr are 
assumed. Here if orm = 5.5 Gyr is adopted according to 
the average formation time obtained by fitting the age- 
redshift relation for four velocity dispersion bins (see Sec- 
tion 2]for details). However, Thomas et al. (2005) found 
that vigorous star formation episodes in massive galaxies 
ofter occur at z ~ 2 — 5, corresponding to a cosmic age 
of less than 3.3 Gyr. The tf orm obtained from the age- 
redshift relations for the quiescent LRGs in this paper 
is larger than that obtained by Thomas et al. (2005), 
which needs further investigation. As the galaxies in 
our sample are all LRGs, the age of the newly formed 
stellar population by major mergers, if significant, is at 
least 1 Gyr, which corresponds to 5 Z ~ 0.1, as the mi- 
gration of galaxies driven by mergers f rom blue cloud to 
red s equence may last ~ 1 Gyr (see iSchawinski et al.l 
I2010D . If a galaxy at redshift z is the remnant of a 
major merger, the age of the younger stellar population 
in it is approximately l/2[iy(z) — t\j(z ~ 0.4)] + lGyr. 
Then a "forged" combined spectrum of galaxies at red- 
shift bin z can be approximately generated by two stellar 
populations, i.e., an old stellar population with the age 
t\j(z) — iform, and a young stellar population with age 
~ l/2[£u(z)-£u (z = 0.4)] + lGyr. The fraction of major 
mergers is 

iHopkins et al. I (|2010T l have shown that the fraction of 
the young population generated by major mergers is 
~ 0.05 - 0.10 in galaxies with > lO n M (see their Fig- 
ure 14). Here we assume it is ~ 0.07. Note that this 
fraction may be an upper limit as the LRGs studied here 
are quiescent ones, in which the star formation gener- 
ated by major mergers might be even less. Therefore, 
the fraction of the young population contributing to the 
combined spectrum is ~ 0.07 x N mrg (z), and the fraction 
of the old population is 1 — 0.07 x N mTg (z). We fit these 
spectra obtained for each redshift by the same method 
as that in Section [3.3.1l and obtain the age for these com- 
bined spectra. Similar to that in Section 21 we obtain the 
best fit of H = 86T^km s" 1 Mpc" 1 by fixing ft m = 0.27. 

From the above calculation, we conclude that Hq may 
be systematically overestimated by up to ~ 20% if us- 
ing the age-redshift relation obtained from the combined 



spectra. Considering this systematic bias, Hq estimated 
from the age-redshift relation in this paper is consis- 
tent with those estimated by other techniques. Note 
also that the systematic bias introduced to the Hq es- 
timate appears not significant for the two sub-samples 
with the highest velocity dispersions, which might mean 
that those very massive LRGs are really quiescent and 
have approximately zero star formation at redshift z < 
0.4. [lb ieiro fc Percivall (|2010l) have pointed out that the 
brightest galaxies show the smallest departure from pure 
passive evolution. Therefore, the most massive LRGs 
with velocity dispersion > 300km s _1 may be efficient 
tools to constrain the cosmological parameters, such as 
Hq, through the age-redshift relations extracted from 
their spectra. 

We also check whether the g-r color of the galaxies 
having the combined spectra in each redshift bin is con- 
sistent with that of the model galaxies that have the same 
spectrum as the combined spectrum in the highest red- 
shift bin. To do this, we extract a model spectrum from 
the GS model grid, whose age and metallicity are the 
best fitting results with ULySS for rcdshift bin 0.36-0.39, 
then we let this spectrum evolve toward the low redshift. 
That is, if the age and metallicity of the combined spec- 
trum are t and Zfc/h a t z=0.375, then at a low redshift 
z, we extract a spectrum from model grid whose age is 
t + St and let the metallicity is fixed to Z-p e /u, where 
St = t\j(z) — t\j(z = 0.375). We repeat this process till 
z=0.045, and denote the spectra obtained from the above 
processes as the "forged spectra". For each combined 
spectrum there is also a model spectrum to fit it. We de- 
note this spectrum as the "model spectrum" . Then we 
calculate the g-r color with SDSS filter of those "forged 
spectra" and "model spectra" with IRAF task SYNPH0T. 
We find that the g-r color of the "model spectra" are 
almost the same as that of the "forged spectra", which 
means that the merger effect is not significant for our 
sample. Furthermore, we pick out those spectrum whose 
S/N > 30 from our sample as a sub-sample to derive 
their physical property with full spectrum fitting. This 
sub-sample has 1386 galaxies, and these are almost all 
concentrated in the redshift range from 0.02 to 0.2. We 
fit these galaxies with SSPs. The initial settings are the 
same as that in Section f3 . 3 . 1 1 Similarly, we extract the 
corresponding spectrum from the GS model according 
to the best fit. Then we calculate the g-r color and find 
that dispersion of the g-r color is very small compared 
with that of the combined spectrum, which means the 
results of the single spectrum are consistent with that of 
the combined spectrum. 

5.4. Model dependence 

In this paper, we use the GalexEV/SteLib model to 
fit the combined spectra, as its wavelength coverage 
is wider compared with the other models provided by 
the ULySS package. For completeness, we also test 
the Pegase-HR/Elodic3.1 model and the Vazdckis/Milcs 
model, and find there is some model dependence. The 
model dependence may be caused by the limitation of 
wavelength coverage as the Pegasc-HR/Elodic3.1 model 
covers the wavelength from 3900A— 6800A, and for the 
Vazdekis/Miles model it is 3540A-7409A. This re sult 
is consistent with that of IVcrkhoda nov et al.l ([2005D . in 
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which they analyzed the photometric data of a large 
sample of elliptical galaxies and found that Hq = 
72± 10km s" 1 Mpc" 1 and 53±10 km s" 1 Mpc^ 1 for two 
different stellar population synthesis models, i.e., PE- 
GASE and GISSEL , respectively. The main ingredients 
of those models are the stellar evolution tracks, the stellar 
library, the IMF, the grids of ages and metallicities, and 
the SFHs, etc. Each model may have different settings 
in one or more of those ingredients. As demonstrated by 
iChen et al.l (|2010l ). at present there are still some differ- 
ences in the output between these different models. The 
model dependence may be avoidable if carefully choosing 
a compatible stellar population synthesis model for the 
problem to be studied. We note here that is also interest- 
ing to test the model given by iMaraston fc Stromback I 
(|2011l) . which is based on the fuel consumption theorem 
and quite different from the other models. However, it is 
not easy to test this model because it is not yet included 
in the ULySS code. 

From the above tests, we conclude that the age- 
rcdshift relation obtained from the combined spectra of a 
large sample of quiescent LRGs can be used to constrain 
the Hubble parameter Hq (and possibly other cosmolog- 
ical parameters if combining with other data sets). If 
a large sample of very massive quiescent LRGs can be 
obtained by the future surveys, such as the Baryon Os- 
cillation Spectroscopic Survey (BOSS), which is less af- 
fected by the systematic bias due to new star formation 
at low redshift, the Hq may be able to be determined 
with substantial accuracy. 

6. CONCLUSIONS 

In this paper, we selected 23, 883 quiescent LRGs from 
the SDSS DR7 in the redshift range from 0.03 to 0.39, 
by setting a threshold of zero emission (at a 2 — a- 
level) of the H Q and [On] lines directly obtained from 
the MPA/JHU catalogue. The quiescent LRG sample 
is divided into four sub-samples according to galaxy ve- 
locity dispersions. For each sub-sample, the spectra of 
galaxies in each of the 12 redshift bins (from z = 0.03 to 
0.39 with a step of 5z = 0.03) are combined together to 
obtain a high S/N combined spectrum. Using the full 
spectrum fitting method, the luminosity-weighted physi- 
cal properties, such as the velocity dispersion, the metal- 
licity and the age, of those quiescent LRGs are obtained 
from the combined spectra by adopting a single popu- 
lation synthesis model, i.e., the GalexEV/SteLib model. 
Using Monte-Carlo simulations, we find that the model 
results arc robust and reliable. We argue that the age- 
redshift relation estimated from the LRG sample could 
be systematic biased because of the contamination from 
a possible younger stellar population formed at z < 0.4 
as consequence of major mergers. This bias is most sig- 
nificant for LRGs with smaller velocity dispersions but 
insignificant for the most massive LRGs. Considering of 
this systematic bias, the age-redshift relation obtained 
from the model fittings is fully consistent with the ex- 
pectations from the ACDM cosmology. 

The Hubble parameter Hq is first estimated by 
using the age-redshift relation obtained from each sub- 
sample, and its value ranges from 89lgkm s" 1 Mpc -1 , 
83^km s- 1 Mpc" 1 , 72+?km s" 1 Mpc" 1 , to 
65^3km s _1 Mpc -1 for the four sub-samples with 



velocity dispersions from low to high, respectively. The 
large value of the Hq estimated from the sub-samples 
with low velocity dispersion is probably due to the 
systematic bias, which can be as high as ~ 20%. 
Using the age-redshift relations obtained from the 
two sub-samples with high velocity dispersions or the 
sub-sample with the largest velocity dispersion, we find 
H Q = 74+|km s" 1 Mpc" 1 or Hq = 65±^km s" 1 Mpc" 1 
if assuming a spatially flat ACDM cosmology, which 
may be less affected by the systematic bias, close to the 
true Hq, and are well consistent with the best estimates 
through other techniques. However, it needs further 
test on whether those most massive galaxies arc truly 
passively evolving or not. 

In summary, we have demonstrated that the age- 
redshift relation of quiescent galaxies can be reliably es- 
timated by using the full spectral fitting method if the 
S/N of their spectra are sufficiently high. We conclude 
that some cosmologies! parameters, such as the Hub- 
ble parameter, can be constrained with considerable ac- 
curacy through the age-redshift relation obtained from 
those most massive LRGs, which is totally independent 
of other methods. With future surveys like BOSS, the 
Hubble parameter may be tightly constrained by the age- 
redshift relation obtained from the most massive quies- 
cent LRGs. 
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